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ABSTRACT 


1/ 

The  present  investigation  expounds  on  the  effect  of  a  ne3r-surface 
cavity,  when  the  solid  surface  is  subjected  to  the  Coulomb  frictional 
loading  of  an  asperity  moving  at  moderately  high-speed.  The  medium  under 
consideration  is  represented  by  a  solid  half  space  which  is  coated  with  a 
thin  layer  of  solid  wear  coating.  The  cavity  in  the  present  report  is 
rectangular.  The  temperature  field  and  its  gradient  in  the  vicinity  of  the 
cavity  result  from  the  traverse  of  the  asperity  over  the  wear  surface  near 
the  cavity.  The  cavity  defect  results  in  a  material  nonuniformity  mathe¬ 
matically  modelled  in  terms  of  the  material  coordinates.  The  resulting 
governing  differential  equation  is  time-explicit  and  transient.  A  general 
finite  difference  formulation  is  developed,  from  which  numerical  solutions 
were  obtained  for  problems  with  a  cavity  at  various  positions  relative  to 
the  surface-layer/ substrate  interface.  Because  of  the  poor  heat  transfer 
characteristics  of  the  cavity,  the  temperatures  in  the  surface  layer  above 
it  are  higher  than  those  in  the  surrounding  region.  This  phenomenon  causes 
3  higher  temperature  gradient,  especially  at  the  trailing  corner  of  the 
rectangular  cavitv. ^There ,  the  direction  of  the  maximum  temperature 
gradient  is  at  a  significant  oblique  angle  to  the  wear  surface.  The 
combined  effect  of  a  higher  temperature  gradient  and  an  oblique  direction 
would  result  in  a  much  larger  shear  stress  in  the  surface-layer/ substrate 
interface.  This  phenomenon  could  lead  toward  delamination  of  the  coating 


near  the  cavitv. 
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Dimensionless  heat  flux  in  ^-direction. 

Dimensionless  heat  flux  in  redirection . 

Dimensionless  temperature  field,  showing  thermal  property  effects. 

Case  2B,  2C,  2D,  2E  -  coated  material  with  a  cavity,  the  top  edge 
of  the  cavity  is  at  the  interface,  showing  thermal 
property  changes  in  comparison  to  Case  4. 

Case  2B  -  The  material  properties  of  the  substrate  are  the  same 
as  the  surface  layer  except  that  k.,  *  2k1 , 
pc0  *  2 pc, . 

Case  2C  -  The  material  properties  of  the  substrata  are  the  same 
as  the  surface  layer  except  that  k0  *  0.5k^, 
pc0  ■  0. Soc^ . 

Case  2D  -  The  material  properties  of  the  substrate  are  the  same 
as  the  surface  layer  except  that  *  k, , 

PC2  =«  0 . 5pCj  . 

Case  2E  -  The  material  properties  of  the  substrate  are  the  same 
as  the  surface  layer  except  that  k0  ■  2k^ , 

pc9  •  PCj . 

Case  4  -  single  material  with  a  cavity,  the  material  is  silicon 

carbide,  L  *  0.04a. 

Dimensionless  heat  flux  in  ^-direction;  n  *  0. 

Dimensionless  heat  flux  in  ^-direction ,  n  ■  0.0^. 

Dimensionless  temperature  field  for  single  material  without  cavity 

and  coated  material  with  cavity,  L  ■  0.04a. 
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NOMENCLATURE 


a  asperity  width 

b  substrate  thickness 

C  specific  heat 

c(t)  distance  from  x^  origion  to  the  center  of  the  asperity 

D  dimensionless  coating  thickness  (H/a) 

d  half  width  of  the  cavity 

e  depth  of  the  cavity 

H  coating  thickness 

k,,k0  thermal  conductivity  of  surface  layer  and  substrate 

4.  — 

L  ligament  thickness 

P(i,j)  center  of  the  finite  difference  cell 

M,N,R,S  surrounding  points  of  P 

P'(x^,t)  pressure  over  the  contact  area 

q(x^,t)  heat  flux  through  the  contact  area 

qQ  average  heat  flux  through  the  contact  area 

0c_^p  heat  flux  from  point  S  to  point  P 

Q*(-q/q^)  dimensionless  heat  flux  through  the  contact  area 

Rj,R0  Peclet  numbers  of  surface  layer  and  substrate  ( R^  *  va Ik 

I  temperature 

U  internal  energy 

v  traverse  speed  of  asperity 

w  width  of  the  total  rectangular  finite  difference  region 

w'  distance  from  x^  origin  to  the  center  of  the  cavity 

£  dimensionless  coordinate  (»x^/a) 

dimensionless  coordinate  (-x^/a) 
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CSTT 


-  v  -  v  * 


C.n 


r i ,IC? 
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coordinates  in  transformed  plane 

thermal  diffusivitv  of  surface  layer  and  substrate 
Coulomb  coefficient  of  friction 

9 

thermal  diffusivitv  ratio  (a"  ■ 

thermal  conductivity  ratio  (B  -  k0/k^) 

dimensionless  temperature  (6  »  Tk, /q  a) 

1  o 

mass  density 

dimensionless  time  (t  =»  vt/a) 
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CHAPTER  I 


INTRODUCTION 

The  present:  investigation  constitutes  a  portion  of  the  general  problem 
of  the  thermo-mechanical  effect  of  a  moving  asperity  over  a  solid  medium 
with  a  solid  coating  layer.  The  asperity,  generally  on  the  order  of  a  milli¬ 
meter,  traverses  over  the  solid  surface  at  a  uniform  moderateiy-high  speed, 
approximately  15  ms  .  and  excites  the  solid  surface  mechanically  and 
thermally  through  Coulomb  friction.  Vhen  there  is  no  flaw  in  the  uniformly 
coated  region,  the  problem  can  be  formulated  with  a  asperitv-fixed  coordi¬ 
nate  system  for  steady-state  solutions.  Ju  and  Chen  [ 1 j  first  solved  for 
the  case  of  a  moderately  thick  coating.  Later  Ju  and  Liu  [2j  extended  the 
general  formulation  of  [1]  to  study  the  thickness  effect  of  the  coating 
layer  for  various  mechanical  and  thermal  impedance  mismatching  between  the 
surface  coating  layer  and  the  substrate.  In  both  cases,  the  results  of 
Huang  and  Ju  [ 3 j  were  verified  that,  at  such  a  moderately  high  speed,  the 
stresses  from  the  thermal  effect  of  the  asperity  friction  are  an  order  of 
magnitude  larger  than  those  from  its  mechanical  traction  effect.  The 
present  investigation,  therefore,  will  concentrate  on  the  temperature  field 
and  the  temperature  gradient  surrounding  a  flaw  in  the  form  of  a  rectangular 
cavity.  This  cavity  is  located  in  the  surface  layer  or  at  the  coating- 
layer  /  substrate  interface. 

The  physical  model  of  asperity  friction  was  based  on  the  experimental 
studies  in  works  by  Archard  [4]  and  Bannerjee  and  Burton  [ 3 ] -  Because  of 
the  three-dimensional  aspect  of  those  observed  "hot  spots",  Huang  and  Ju  [3j 
formulated  the  problem  in  three-dimensional  theory  of  the  thermo-mechanical 
effect  owing  to  the  Coulomb  friction  of  moving  asperities  over  a  single 


solid  wear  medium.  In  its  application  to  the  problem  of  a  sir.  .la  high  speed 
asperity,  it  was  observed  that:  (i)  at  an  asperity  traversing  speed  of 
10-15  ms  *,  the  thermal  stress  as  a  result  of  high  temperature  gradients  is 
eight  times  the  stress  due  to  mechanical  loading,  (ii)  at  such  asperity 
speed,  the  maximum  tensile  stress  is  the  largest  at  a  depth  of  about  one 
tenth  the  asperity  size,  i.e.,  at  a  depth  of  the  order  of  100  pm,  and  (  iii) 
while  the  surface  temperature  is  the  same  for  the  three-dimensional  model  as 
compared  with  that  of  two-dimensional  model  16',  the  resulting  maximum 
stress  is  six  times  higher  for  the  three-dimensional  model  due  to  much 
higher  temperature  gradient.  Therefore,  for  the  thermo-mechanical  analysis 
of  the  effect  of  an  imbedded  flaw,  it  is  primarily  important  to  obtain  first 
a  temperature  field  solution. 

In  the  previous  works  of  the  moving  asperity  problem,  the  analyses 
dealt  with  basically  uniform  solid  medium,  that  is,  the  material  and 
asperity  properties  are  invariant  in  the  direction  of  the  asperity  motion. 

In  such  a  case,  the  time  effect  can  be  rendered  implicit  in  the  Fourier  and 
'Javier  equations  bv  using  a  coordinate  system  fixed  to  the  asperity  (called 
the  convective  coordinate  system)  traversing  therewith  [i,2,3i.  The 
resulting  steady-state  solutions  are  not  explicitly  time-dependent  in  such 
coordinates.  When  a  flaw,  such  as  a  cavity,  exists  in  the  medium,  unifor¬ 
mity  in  the  direction  of  asperity  motion  no  longer  exists.  A  coordinate 
system  fixed  to  the  cavity,  or  to  the  material,  (referred  to  as  the  material 
coordinates)  must  be  employed  by  necessity.  As  the  asperity  traverses  over 
the  surface  under  which  the  cavity  lies,  the  induced  temperature  field  has 
to  be  an  explicit  function  of  time.  The  solutions  are  transient.  The 
steadv-state  solution  of  the  temperature  field  resulting  from  the  moving 


asperity  over  an  uniform  medium  has  been  dealt  with  bv  several 


ors 


1.1,2,3,6,7,8,9,10],  whereas  the  transient  solution 
lacking . 
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for  a  flawed  medium  is  •[-! 
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CHAPTER  II 


ANALYTICAL  MODEL 


WTiile  the  three-dimensional  model  gives  a  more  accurate  description,  it 
is  considered  that,  for  the  purpose  of  studying  the  effect  of  a  cavity,  a 
two-dimensional  model  is  adequate,  with  considerably  simpler  mathematical 
operation.  The  general  geometry  is  shown  in  Figure  1. 


2.1  Governing  Mathematical  Equations 

The  mathematical  model  is  a  semi-infinite  body  with  a  thin  coating 
layer  and  a  rectangular  cavity  in  the  neighborhood  of  the  coating- 
layer  /  substrate  interface.  The  half  space  surface  is  subjected  to  the 
frictional  heating  of  a  moving  asperity  over  the  surface. 

The  governing  equations  are  as  follows:  In  the  surface  layer, 

0  <  n  <  D,  denoted  by  the  superscript  (1), 


,2  (1)  ,  2  ( 1 )  , 

3  o  3d  „  3  o 

-  +  - -  „  Ri  _ 


(  1) 


(1) 


35“  3rf  *  3t 

In  the  substrate  region,  D  <  n  <  °°,  denoted  by  the  superscript  (.2), 

(2) 


32o(2)  32o(2) 


3b' 


(2; 


3C 


3n‘ 


3t 


where  o(“Tk../q  a)  is  the  dimensionless  temperature,  (  5 ,  n)  (  *x . /a)  are  the 
lo  1 


dimensionless  coordinates  in  the  direction  opposed  to  the  asperity  motion, 
and  the  depth,  direction  respectively,  t(avt/a)  is  the  dimensionless  time. 


D  =  H/a  the  dimensionless  coating  thickness,  *  va/r,  the  Peclet  numbers 


in  the  coating  layer  (i»l)  and  in  the  substrate  (i»2),  and  Q*(aq/qo)  is  the 


dimensionless  heat  flux  from  the  asperity  friction  normalized  by  the  average 


. 1  % 
.'•'.‘AV 


heat  flux  qQ  in  the  contact  area,  v  is  the  asperity  traversing  sree d,  a  is 
the  asperity  width  in  the  traversing  direction,  ir  is  the  thermal  diffusiv 
ity  in  the  region  (i). 

The  temperature  field  6  must  satisfy  the  following  conditions: 

1'"  the  initial  condition, 

O(i)(£,n,0)  =0,  i  *  1,2  (3) 

2''  the  boundary  conditions: 

•>  o 

(i)  at  infinity,  £”  +  n”  -*  the  regularity  conditions  are 

o(1),o(2)  -  0, 

<ii)  in  the  asperity  contact  surface,  c(t)/a  <  ^  <  (c(t)/a)+l,  n  *  0,  the 
heat  input  through  the  boundary  is  the  rate  of  the  frictional  energy 


"  —  *  °*  =  *'%' 
dr)  o 

where  q  *  )ifvP'  ,  u.P'  is  the  asperitv  friction  force  in  the  contact  regi< 
r  r 

(iii)  outside  the  contact  surface,  5  <  c(t)/a  or  C  >  (c(t)/a)+l,  n  *  0, 
without  loss  of  generality,  the  surface  is  adiabatic 


3“  at  the  layer/ substrate  interface,  n  *  D,  the  continuity  conditions  for 
the  heat  flux  must  be  satisfied. 

o(1)  -  *(2) 


3n 


where  8  (*k0/k^)  is  the  thermal  conductivity  ratio, 
'■*'  adiabatic  conditions  at  cavity  boundaries. 


—  ■  0,  at  £  -  id/a,  L  £  n  £  (L+e)/a, 


ASPERITY 


.(2),  .  .  , 
<b  ( l,  j  ,n) 


$ 


a 


1 

*$1 


1 


a2r^  1 4>(  2  *  (  i-1  ,;j  ,n-l !  +  b*'  ^  *  (  i+1 ,  i  .n-1 )  J  + 


?  (?) 

+  [1  -  2a“(  r1  +r9 )  J<6  “  (i,j,n-l)  + 


2  (  ? )  ( 2 ) 

+  a  r^[6  “  (i,j-l,n-l)  +  o  ( i, j+1 ,n-l ) J . 


r  f  * 

where  r^  ■  At/ ( R^ • AE") ,  r0  *  At/ ( R^ 'An- ) ,  c<“  ■  r^/ir^,  and  (i,j,n)  denotes 


the  two  spatial  indices  and  the  time  step  respectively.  The  stability- 


criteria  for  the  explicit  scheme  can  be  shown  to  be 


1  —  2 ( r j  +r0  )  2  0 , 


1  -  2a“( r^+r., )  £  0. 


If  ct“  is  less  than  1,  Equation  (12)  is  the  stability  criterion,  otherwise 


Equation  (13)  is  the  stability  criterion  for  the  explicit  scheme. 


Based  on  the  previous  works  l 1 , 2 , 3, 6 , 10 J ,  we  know  that  the  high 


temperature  and  the  high  thermal  stresses  occur  in  the  region  near  the 


asperity.  Therefore,  in  that  region  one  must  use  very  fine  mesh  to  obtain 


good  results.  But  for  those  locations  far  away  from  the  asperity,  one  can 


use  coarse  mesh  in  order  to  save  computing  time.  This  non-uniform  mesh  can 


be  taken  care  of  by  applying  the  general  coordinates  transformation  [131 


(Fig.  3).  The  derivatives  of  a  function  f(E,n>  in  the  transformed  plane 


(  E .  n )  are : 


f .  -  ( n-  f T  -  HT  f-)  /J  , 
E  n  E  E  n 


f  -  (E7  f-  -  E-  fT)/J  , 
n  E  n  n  E 


f__  -  (  o-2  f  r r  ~  2rrr  n-  fp-  +  f - )/  j" 

EE  n  EE  E  n  En  E  nn 


+  Un;2  npp  ~  -h7  n-  np-  +  ivp2  n--)  (Z~  £j  -  Zj  £-) 
n  EE  E  n  En  E  nn  n  E  E  n 


+  <n-2  Epp  -  2n7  n-  Ep-  +  n72  E-->  (nrf;-  n-  f7>]/J3, 
n  EE  E  n  En  E  nn  E  n  n  E 


»  *  '  M  *  •  «  '  '  -  *  * 
■  .  1 

:  *  -  ■*.  -  *.  -  -  -  *  .  <- 


v  -  .  •*  .*• 


f  «  tc-  t--  -  2C7  e-  f?-  +  £7' 

nn  n  ££  £  n  £n  £  nn 

+  l(£-“  t)-?  -  2£-  £—  n?-  +  £p2  n--)  <  £-  fr  -  £7  f-) 

n  ££  £  n  £n  £  nn  n  £  £  n 

♦  <£-2  £?r  -  ^£7  £-  £r~  +  ^  C~  >  <  IT  f~  -  H~  f  p  >  I  /  J 3 .  (  17* 

n  ££  £  n  £n  £  nn  £  n  n  £ 

where  J  ■  £7  p—  -  hT  is  the  Jacobian  of  the  transformation,  the 

subscripts  (£,n,£,n)  denotes  partial  derivatives  in  those  coordinates 

respectively. 

The  heat  conduction  Equations  (1,2)  in  the  transformed  plane  can  be 
written  as 

,,  (i)  »,  (i)  ,  .(i)  ,  (i)  ,  (i).,.2  „  ii) 

(A.  0CC  -  2A-  br  +  A-  0  +  A,  0  +  A-  0_  )/J  ■  R.  0 

1  ££  2  £n  3  nn  4  n  0  £  1 

i  -  1,2  (13) 


where 


0  0 

A  -  £-'  +  n-“  , 

inn 


A:  '  5e  *  n5  %  ' 


0  9 

A3  ‘  *  V  • 


-  1  *5  A6 

-  £7  A. ) / J  , 

£  > 

(  22) 

-  '  ^  A7 

-  n-  A  ) / J  , 
n  0 

(  23) 

1  He 

-  2a2  £^~  +  a3 

^  • 

(  24) 

*  A1  n?f 

-  -A2  n£n  +  A3 

nnn 

(25) 

At  the  outer  boundaries  of  the  rectangular  region  (excluding  the 
surface),  0  *  0  is  the  nominal  value.  The  remaining  conditions,  on  the 
surface,  the  cavity  boundaries  and  at  the  interface,  shall  be  incorporated 
with  an  energy  balance  scheme  114). 


2.3  Energy  Balance 

The  cavity  boundaries,  the  moving  asperity  and  the  interface  of  the 
medium  are  taken  care  of  with  the  use  of  the  energy  balance  method. 

( i)  Energy  balance  at  the  interface  (see  Fig.  4). 

For  the  material  1  (surface  layer),  the  heat  fluxes  toward  the  point  P 

at  the  interface  from  material  points  M,  R  and  S  in  the  surface  layer  are 


.„?  ■  r  1-1  Ill~1> 


mw:,  T,i*l-1’'-  Tlt't' 

i  •  Ax 


QS-P  *  kl(^> 


T  (  i  .  i  - 1  )  -  Ki.j) 


(27) 


(23) 


For  the  material  2  (substrate),  the  heat  fluxes  toward  the  interface 
point  P  from  material  points  M.  R  and  N’  in  the  substrate  region  are 


0  *  k  (Av/'))  ~  > — 71— II iaJ-l 

■M-P  2  •  Ax 


0  *  k  (  Av / 0  )  ^  J  }  ~ 

K2  '  Ax 


■R-*P 


T(  i.  :  +  l  )  -  T'  i,  j  ) 


Q,.  n  *  k-( Ax) 

N-*P  2  Ay 


(29) 

(  30) 

(  31 ) 


where  Q  is  the  heat  flux,  indexed  by  the  flow  direction. 

The  total  heat  flux  goes  in  the  interface  point  P(i,j)  is 

Q  •  tic  +k,>  (Av/2)  fli  1~l  ?  ~  +  T(l*1,J>  ~ 

xsum  1  2  L  Ax  Ax  J 


+ 

(32) 


+  k  (Ax)  T(  +  k  (Ax) 

1  Ay  2 


T( i . ]+l )  -  T( i .  i ) 
Ay 


The  rate  of  change  of  internal  energy  U  in  the  time  interval  At  and 


in  the  neighborhood  of  P(i,j)  is 


/  / 


For  conservation  of  energy,  U  ■  Q  ,  one  obtains  the  equation  for 
the  continuity  condition  at  the  interface  point  P(i,j)  at  the  time  step  n 
is 


T( i, j ,n) 


!2(k1+k2)  At  2(k1+k2)  At  ' 

*  “  ‘  2  ~  *  ~ 

( p^c^p^c.,)  Ax  ( p^c^+p9c9 )  Ay"_ 


T(i,j,n-1)  + 


<  kj+k,) 


At 

— -  [Tli-l.j.n-l)  +  T( i+1 , j ,n-l ) J  + 


( p, c. +p0c9  )  Ax* 


11 
2k 


At 


( p^c1 +p9c9 )  Av' 


•  — ~  T( i, j-1 ,n-l )  + 


(p^c^+p^c,,)  Ay' 


At 

— ^  T< i,  j+1 ,n-l )  . 


( 34) 


Equation  (34)  in  dimensionless  form  is 


0( i, j ,n) 


[“ 


2< 1+8) 


(1  +  B/oT) 
(  1  +  3  ) 


1  rl 


2(  1+3) 


2.  *2] 


0( i, j , n-1 ) 


( 1 +3 /a 
o 


( 1+3/cO 

T~  rjiO(  i-1 , j ,n-l )  +  d(  i+1 , j ,n-l) ] 

2B 


(  35) 


->  ‘•9 

( l  +  B/a“ )  ' 


r„o( i , j-1 ,n-l ) 


(  1+3/cT) 


r9o(  i, j  +  1 ,n-l )  , 


where  B  3  k0/k1 . 

(ii)  Energy  balance  on  the  cavity  boundaries  (Fig.  5).  The  heat  fluxes 
toward  the  boundary  point  P(i,j)  from  the  material  points  M,  S  and  N  are 


^M~*P 

Qs-*P 

Q.V-*P 


k(Av)  — ■;L~1  ’■! }  T  T(  1i.1  ’ 

Ax 

k(  Ax/2  )  Tf  ~  Tf  i’'l)  , 

Ay 

k( Ax/2)  — +  I  1  "  T( i'^ -  . 


(  36) 

(37) 

(  38) 
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Figure  5.  Energy  balance  on  the  cavity  boundari 


The  total  heat  input  at  P(i,j)  is  then: 

I.*.  ■  Vp  *  <3~P  *  Vp  '  ^  T<1~1'1'ak  111,11  + 


+  k( Ax/2)  -(  1:Jtl }  ~  T(  +  k(Ax/2> 

Av 


T(  i. j+1  )  -  T< i.  i  ) 
Av 


The  rate  of  internal  energy  increase  in  the  neighborhood  of  P(i,j)  is 


U  *  oc  (0.5  Ax  Av) 


T( i. j .n)  -  T( i. i .n-1 ) 


From  conservation  of  energy,  one  can  obtain  the  temperature  at  the  boundary 
point  P(i,j)  at  the  time  step  n: 


A  t 

i,j,n)  =  T(i.j,n-1)  +  - 

oc 


T(i-l,j.n-l)  -  Tti.j.n- 


♦  0>5  Kid-U-D  - 


2T(  i  ,  i  ■  n  - 1 )  +  T(  i.  j-tl  .n-1  )  j 
Av2 


Equation  (  H  )  in  dimensionless  form  is 

o(i,j,n)  *  o(i,j,n-l)  +  2r, C o ( i-1 , j ,n-l )  -  o(i,j,n-l)J  + 

+  r 0 C o( i, j-1 ,n-l )  -  2o( i , j , n-1 )  +  0( i  ,  j  +  1 ,n-l ) J  .  (12) 

Similarly,  other  points  on  the  cavity  boundaries  can  be  obtained. 

(iii)  Energy  balance  at  the  corner  of  the  cavity. 

The  points  at  the  four  corners  of  the  cavity  are  singularities  because 
at  each  of  those  four  points  there  are  two  boundary  conditions, 

<),T  ■  3^1  *  0,  with  only  one  unknown  T.  However  by  applying  an  energy 
balance  scheme,  one  can  resolve  such  problems  at  the  comer  (see  Fig.  6). 
The  heat  fluxes  toward  the  corner  point  P(i,j)  are: 


Qm-p  " 


T( i-1 . i )  -  T< i, i  ) 


QN_p  *  k<  Ax  /  2 )  Tf  1;J- 


w 


*  -  -  r.  . V  .  -u  -VrSv  •  -  _•  j  ‘  •  •  -  '•  *  *  v  *• 
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QR"*P  *  k(Ay/2)  ^  T<1,1~ 


Qs-p  = 


T( i. 1-1  )  -  T( i.  i  ) 


The  total  heat  input  at  P(i,j)  is  therefore 


Q  *  k( Av ) 
sum 


V.'r'-H.  -  +  k!Ax/2)  ILL, 


t! )  -  Tt i,  i  ) 
Av 


,  „  T(  i+1 .  i  !  -  ■r(  i .  \  )  T!  i  '■ 

-  k!  Av  /  2 )  — ,  — - : - •-— -—  *  k'A:;> 


1  )  -  I(i,' 


The  rate  of  increase  in  internal  er.errv  ia 


U  =  o c <  0.5  Ax  Av  +  0.25  Ax  Av ) 


T  (  i ,  j  .  n  >  - 


From  conservation  of  energy,  the  equation  for  boundary  con; 
the  corner  point  Pti.j)  is  therefore 


Tti.j.n)  =  T.'i.j.n-l)  +  (A/3)  (k/oc)  (At/Ax")  (T(i-l,n-l) 


-  (3/2)  Ki.j.n-l)  +  (1/2)  T(  i+1 ,  j  ,n-l )  ]  r 


r  (A/3)  (k/'pc)  (At/Av~)  { Tt  i  .  j -1 ,  n-i  )  - 


(3/2)  Tt i  ,  j , n-1 )  +  (1/2)  Tt i , j +1 , n-1 ) ] 


location  ' -9)  in  dimensionless  form  is 


o(i,j,n)  ®  o(i,j,n-l)  +  <  A/ 3 )  r  ^  [  <M  i-1 ,  j  ,  n-1  )  -  1.5  o(i,j,n-l)  + 


+  0.5  0( i+1 , j ,n-l ) ]  +  ( 4/ 3 ) r0 [ <M  i , j-1 ,n-l )  - 


-  1.5  o(i,j.n-l)  +•  0.5  o(  i  ,  j  +  1 ,  n-1 )  ]  . 


ST 


Similarly,  the  boundary  conditions  at  the  other  corner  points  of  the  cavity 
are  obtained. 


V.' 


*  . 

(•Y, 


(iv)  Energy  balance  on  the  surface  boundary  (under  the  moving  asperity  at 
the  time  step  n) . 

For  the  explicit  scheme,  the  time  step  is  limited  by  the  stability 
criterion.  As  a  result,  the  moving  asperity  at  some  time  may  not  be  right 
above  the  grid  points.  To  alleviate  this  situation,  one  can  also  use  the 
energy  balance  method  to  describe  the  boundary  condition  (see  Fig.  7).  The 
heat  fluxes  toward  P  from  material  points  M,  N  and  R  are 


q_m_p  ■  uuy/2) 


(  51  ) 


Vp  *  k(Ax) 


T(  i , 2  )  -  T(  i. 1 ) 
Ay 


Qr-,o  *  k(Ay/2)  — i--U-1-)-"-Tj-i-"U 


(  52) 


(53) 


The  exterior  heat  flux  in  the  neighborhood  surface  of  the  boundary  point 
P(i,j),  which  is  under  the  asperity  at  the  time  n.  is 


Q  _  =>  q  (  0 . 5  Ax  +  h)/unit  thickness 
ext-*P 


where  h  is  less  than  Ax/2.  The  formulation  thus  takes  care  of  all  cases 
when  the  asperity  end  points  do  not  fall  on  a  material  grid  point  on  the 
boundary.  The  total  heat  input  at  P  is  therefore 


sum 


k  fUv/2,  m-lalLz  Tf  i  >— ■  +  (Ax)  T(1> 

L  Ax  Ay 


1  ) 


.  ,  ,  T(i+l.l)  -  Tf i , 1 ) 1  .  ,  .  .  ,  , 

+  (Ay/ 2)  - — - 1 —  +  q(0.5  Ax  +  h)  . 


(  55) 
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The  rate  of  increase  of  the  internal  energy  in  the  neighborhoc 

-  .  .  .  T(i.i.n)  -  T(i.i.n-l) 

l  *  pc  (0.5  Ax  Av)  - J ^ -  . 

At 

The  condition  on  the  surface  boundary  with  the  use  of  the 
of  energy  becomes 


■i  of  P(i,j)  is 
(  56) 

conservation 


T(i.l.n)  -  T(i.l.n-l)  +  (k/pc)  (At/Ax2)  [T( i-1 , 1 ,n-l )  -  2T(i,l,n-l)  + 

+  T!  i-cl ,  1 ,  n-1 )  ]  +  (2k/ pc)  (At/Ay2)  [T(i.2,n-l!  -  T(i,l,n-D]  + 
2a<  0. 5  A::  +  h)At 

~  — — — -  .  (  o  /  ) 

pc  Ax  Ay 

Equation  ( 57 )  in  dimensionless  form  is 

o'i.i.n)  *  o(i,l,n-l)  -t-  r,  [o(  i-1 , 1  ,n-l )  -  2o(i,l,n-l)  +  o(  i+1 , 1 ,  n-I  )  ]  + 

x 

+  2r2lo(i,2,n-l>  -  <b(  i ,  1  ,n-l )  ]  +  d+h’)  ,  (53) 

where  h'  ■  h/(Ax/2). 

Equations  (10,11.35,42,50,58)  constitute  the  general  formulation  of  the 


problem  with  a  complete  set  of  difference  equations  for  the  solutions  of  the 
discrete  temperature  field  (d(i,j,n))  at  some  specific  time. 


CHAPTER  III 


NUMERICAL  RESULTS 

Numerical  results  are  obtained  by  using  the  non-uniform  rectangular 
mesh  corresponding  to  different  cases  of  material  properties  and  geometry. 
For  the  surface  layer  of  Silicon  Carbide,  k^  3  1.047  J/cm.“C.s, 

O 

k',  3  0.4?  cm"/s,  and  c ,  3  712  J/kg.“C.  For  the  substrate  of  Aluminum, 

k0  =  2.02  J/cm.^C.s,  3  0.961  cm”/s,  and  c?  =  917  J/kg.JC.  The  other 

__  1 

numerical  parameters  on  the  asperity  and  the  cavity  are:  v  3  15  ms  *, 
v  3  4a,  H  3  1.2a,  £  3  1.9a,  d  3  0.3a,  e  3  0.5a,  a  3  1mm,  the  smallest  A£  and 
An  are  0.02  and  0.01  respectively,  and  At  =  0.01.  At  the  limiting  case  of 
no  cavity,  the  maximum  dimensionless  temperature  at  surface  of  the  coated 
media  was  found  to  be  0.124  by  using  the  transform  method  [lj.  The  result  at 
the  same  point  by  the  current  finite  difference  formulation  is  0.123.  The 
error  is  less  than  l'S>.  The  numerical  scheme  is  therefore  confirmed  by  the 
benchmark  problem. 

The  solutions  for  a  single  material  with  and  without  a  cavity  would 
then  be  compared  with  two  limiting  cases.  For  the  first  case,  the  cavity  is 
located  entirely  in  the  surface  layer.  Figure  8a.  In  the  second  case,  the 
top  edge  of  the  cavity  is  at  the  layer/ substrate  interface.  Figure  8b.  The 
solutions  for  the  single  material  without  and  with  a  cavity  are  designated 
as  the  third  and  the  fourth  cases  respectively,  included  for  the  purpose  of 
comparison.  In  Case  1,  three  different  values  of  the  ligament  thickness: 
0.04  (Case  1A) ,  0.06  (Case  IB)  and  0.1  (Case  IE),  are  used  to  illustrate  the 
effect  of  the  ligament  volume.  The  temperature  fields  at  two  depths,  for 
all  three  Cases  1A,  IB  and  IE  are  shown  in  Figure  9.  In  the  figure,  the 
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cavity  width  is  from  5  *  1.6  to  2.2.  The  asperif  position,  a",  t'.-.e  dimen¬ 
sionless  time  t  “  1.04,  representing  the  worst  case,  is  from  £  =  1.2  to  2. 
The  relative  positions  of  the  asperity  and  the  cavity  is  shown  in  Figure  10. 

The  Case  lA  then  is  compared  with  Case  2  of  the  sane  ligament  thickness 
for  which  the  top  edge  of  the  cavity  is  at  the  layer/ substrate  interface. 

The  effect  of  the  relative  position  of  the  cavity  to  the  interface  is  shown 
in  Figure  11.  It  is  noticed  that  when  the  top  edge  of  the  cavity  is  3t  the 
interface,  the  temperature  field  in  the  region  immediately  on  the  trailing 
edge  of  the  asperity  will  be  affect  by  the  substrate  material. 

The  effect  of  heat  capacity  3nd  thermal  conductivity  of  the  surface 

layer  for  Case  lA  is  shown  in  Figure  12.  The  figure  shows  the  original 

value  as  Case  lA.  Case  1C  represents  a  reduction  of  thermal  capacity  of  the 

surface  layer  by  half.  Case  ID  shows  the  result  of  increase  in  thermal 

conductivity  of  the  surface  layer  by  75”i.  The  thermal  conductivity  of  the 

surface  layer  is  shown  to  have  little  effect  cn  the  nondimensional  surface 

temcerature.  But  the  real  temperature  field.  T  ■  q  ao/k, ,  is  lowered  with 

o  1 

increasing  thermal  conductivity  k, . 

Figures  13  and  11  illustrate  the  effect  of  a  cavity  on  the  direction  of 
heat  flux.  The  figures  show  the  nondimensional  heat  flux  components  in  £ 
and  n  directions  of  single  material  without  cavity  (Case  3)  and  a  layered 
medium  with  a  cavity  (Case  lA) .  From  the  figures,  it  is  observed  that  with 
no  cavity  the  heat  flux  at  £  *  2.2,  and  n  ■  0.0^  has  a  magnitude  of  0.7  at 
an  angle  of  82"'  to  the  wear  surface,  with  cavity,  at  the  same  location,  the 
magnitude  is  increased  to  1.5  at  an  angle  23“’  to  the  wear  surface.  Hence, 
the  existence  of  the  cavity  will  increase  the  heat  flux  tremendously, 
especially  in  the  £  direction  near  the  upper  trailing  comer  of  the  cavity. 
Figures  13  and  14  demonstrate  not  only  an  increase  in  magnitude  of  the  heat 
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flux,  hence  the  temperature  gradient,  but  also  the  flux  at  a  me r r  oblique 
angle  to  the  wear  surface. 

Ju  [ 1  - i  had  studied  the  effect  of  thermal  properties  of  a  single 
material  subjected  to  the  high-speed  asperity  excitation.  It  was  pointed 
out  that  thermal  conductivity  (k)  and  thermal  capacity  (pc'  are  the  para¬ 
meters  controlling  the  temperature  field.  For  layered  media  similar  effect 
was  concluded  bv  Ju  and  Liu  l 2],  For  the  case  of  a  layered  medium  with 
cavirv,  the  thermal  property  variation  in  the  coating  layer  can  be  accord¬ 
ingly  extrapolated.  It  is  the  effect  of  the  substrate  in  the  neighborhood 
of  the  cavity  that  would  be  influential  to  the  temperature  field  in  the 
critical  region.  The  effect  of  thermal  property  variation  for  the  substrate 
is  therefore  studied  numerically  for  the  Case  2,  for  which  the  coating/ 
substrate  interface  is  at  the  top  edge  of  the  cavity.  For  this  case,  the 
thermal  properties  of  the  substrate  will  be  immediately  influential  of  the 
temperature  field  in  the  vicinity  of  the  top  trailing  corner  of  the  cavity. 
Far  the  purpose  demonstration  of  the  individual  effect,  a  benchmark  case  is 

chosen  for  comparison  that  both  the  coating  and  the  substrate  are  of  silicon 

2 

carbide,  ( k,  *  k7  =*  1.047  J/cm.JC.s,  ir  =•  *  0.49  cm“/s, 

3 

pc^  *  pc,,  a  2.137  J / cm  .C),  designated  as  Case  4.  Figure  15  illustrates  the 
temperature  field  near  surface  and  at  the  coating/ substrate  interface  for 
cases  with  marked  changes  in  thermal  properties  from  those  given  in  Case  4. 
Case  2B  shows  no  change  in  the  substrate  diffusivity,  but  both  the  thermal 
conductivity  and  the  thermal  capacity  are  doubled.  The  ensuing  improved 
conductivity  and  capacity  in  the  substrate  allow  a  significant  heat  flaw 
into  the  substrate  thus  a  high  temperature  gradient.  The  Case  2C,  at  the 
same  diffusivity  but  with  both  the  thermal  conductivity  and  the  thermal 
caDacity  halved,  3how3  a  reduced  he3t  flow  into  the  substrate  thus  a  low 


but  from 


temperature  gradient.  Cases  2D  and  2E,  with  doubled  diffusivity. 
half  capacity  and  double  conductivity  respectively,  showed  reduced  and 
increased  heat  flow  into  the  substrate  respectively.  The  heat  flux  being 
proportional  to  the  temperature  gradient,  is  illustrated  in  Figures  16  and 
17  for  the  surface  region  and  at  the  interface.  It  is  to  be  noticed, 
however,  that  the  temperature  gradient  and  the  heat  flux  are  not  of  identity 
especially  across  dissimilar  materials.  The  normal  component,  to  the 
interface,  of  the  heat  flux  must  be  continuous  but  not  the  temperature 
gradient.  The  tangent  component  of  the  temperature  gradient  must  be  con¬ 
tinuous  across  the  coating/substrate  interface.  The  tangent  component  of 
the  heat  flux,  however,  is  generally  not  continuous  across  the  interface. 

In  the  numerical  computation,  using  the  method  of  energy  balance.  Equation 
(32),  the  x^-component  of  heat  flux  at  the  coating/ substrate  interface, 
Figure  17,  is  the  average  of  the  heat  flux  components  in  the  coating  and  the 
substrate,  that  is  -0 . 5( k^+k0 ) 3T/ 3x^ . 

Figure  13  shows  the  transient  temperature  for  Case  1A  (cavity  in  the 

coating  and  ligament  thickness  of  0.04)  in  comparison  to  the  case  of  a 

single  material  without  cavity  (Case  3).  The  dimensionless  temperatures, 

0  =*  Tk, /q  a,  Dlotted  against  the  time,  x  =  vt/a,  at  surface  and  at  the 
I  o 

ligament  depth,  n  *  0.0^,  for  the  position  £  *  2.2  where  the  temperature  is 
maximum  in  the  vicinity  of  the  cavity.  It  is  shown  that  before  the  asperity 
reaches  the  point,  the  temperature  is  low.  Then  the  surface  temperature 
increases  and  reaches  a  maximum  when  the  asperity  just  passes  over  the 
trailing  edge  of  the  cavity.  The  temperature  at  the  trailing  corner  of  the 
cavity  lags  in  time  for  it3  peak  value. 
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CHAPTER  IV 


CONCLUSION 


The  transient  solution  to  a  multiple-boundary  problem  is  always 
complex.  In  the  current  problem,  the  totality  of  the  boundary  is  a  union  of 
the  wear  surface,  the  cavity  boundaries,  the  material  interface  and 
infinity.  When  the  transient  temperature  solution  is  to  be  extended  in  the 
near  future  to  the  transient  thermo-mechanical  stress  solutions,  the  com- 
plexitv  can  only  be  anticipated.  The  numerical  solution  becomes  a  necessity 
for  analyses  in  special  cases.  In  the  present  report,  it  has  been  demon¬ 
strated  that 


1.  the  transient  governing  differential  equations.  Equations  (1,2),  can 
be  formulated  into  difference  equations.  Equations  (10,11); 

2.  variable  grids  (£,n),  made  necessary  due  to  strong  local  effect,  can 
be  transformed  into  a  uniform  grid  space  (^,n); 

3.  boundary  conditions  in  temperature  and/or  heat  flux  can  be  expressed 
through  energy  balance,  thus  avoiding  the  problem  of  singular  points 
such  as  those  of  the  cavity  corners; 

4.  the  numerical  solution  can  be  tested  by  benchmark  with  a  known 
analytical  solution,  showing  an  accuracy  of  within  1 %. 

Like  most  numerical  solutions  of  field  problems,  functional  relationship  is 
not  obtainable  without  volumnous  computations.  However,  significant  con¬ 
clusions  can  be  reached  through  a  few  data  points. 

The  numerical  variation  in  the  thermal  properties  of  the  coating  layer 
confirms  the  effects  of  the  thermal  capacity  (pc)  and  the  thermal  conduc¬ 
tivity  (k),  Figure  12,  as  observed  in  Huang  and  Ju  [10).  Specifically  for 
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Che  ca3a  when  there  is  a  cavity  near  the  wear  surface,  we  can  ;raw  the 
following  conclusions: 

1.  Because  of  the  discontinuity  in  heat  transfer  across  the  void,  the 
temperature  will  rise  higher  in  the  ligament  region  than  the  case 
without  the  cavity,  Figures  9  and  II. 

2.  The  temperature  rise  is  inverse  to  the  ligament  volume,  Figure  9. 

3.  When  the  interface  is  at  the  ligament  depth,  the  thermal  property  of 
the  substrate  will  influence  the  temperature  gradient.  Figure  15. 

4.  Because  of  the  necessary  heat  transfer  in  the  lateral  direction,  the 
heat  flux  will  be  at  a  large  oblique  angle  to  the  wear  surface.  In 
the  case  of  a  layered  medium  without  a  cavity  the  heat  flux  is  in 
the  direction  approximately  9U  from  the  perpendicular  to  the  wear 
surface  at  a  depth  of  0.1a.  If  the  coating/ substrate  interface  is 
in  the  neighborhood,  the  resulting  delaminating  shearing  stress  is 
up  to  30 of  the  maximum  principal  stress,  Ju  and  Liu  (2j.  With  the 
presence  of  a  cavity,  not  only  the  magnitude  of  the  temperature 
gradient  increases,  but  also  the  direction  of  the  maximum  tempera¬ 
ture  gradient  is  rotated  to  a  more  oblique  angle  to  the  wear  sur¬ 
face,  Figures  10  and  11.  Shear  stress  in  the  layer/substrate 
interface  is  expected  to  be  huch  higher.  The  temperature  field 
solution  thus  points  to  the  possibility  of  a  shear  delamination 
between  the  coating  layer  and  the  substrate,  because  of  the  presence 
of  a  cavity  flaw. 
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